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The Advanced LIGO and Advanced Virgo ground-based gravitational-wave (GW) detectors are projected to 
come online 2015—2016, reaching a final sensitivity sufficient to observe dozens of binary neutron star mergers 
per year by 2018. We present a fully-automated, targeted search strategy for prompt gamma-ray counterparts in 
offline Fermi-GBM data. The multi-detector method makes use of a detailed model response of the instrument, 
and benefits from time and sky location information derived from the gravitational-wave signal. 



1. Introduction 



Compact binary coalescence (CBC), such as the 
merger of two neutron stars (NS) or black holes (BH) , 
remains the most highly anticipated gravitational- 
wave signal for ground-based gravitational-wave de- 
tectors. The second-generation Advanced LIGO 
(4km-baseline interferometers in Hanford, WA and 
Livingston, LA) [Harry 2010| and Advanced Virgo 
(3km interferometer in Cascina, Italy) |Accadia et al. 
2012] detectors are projected to come online 2015- 
2016, and reach a final sensitivity sufficient to observe 
dozens of NS-NS mergers per year by 2018. Together 
with the 600m GEO-HF detector iWillke et al.|2006l, 



they will form a word-wide network of gravitational- 
wave interferometric detectors. The network will be 
joined by the Japanese 3km cryogenic KAGRA detec- 
tor Somiya 2012 and a proposed third LIGO-India 



detector around 2018-2020, increasing overall sensi 
tivity and sky-coverage, while also imp roving sky lo 
calization and waveform reconstruction Schutz|2011 



Gravitational-waves have yet to be directly ob- 
served. Our confidence in the existence of compact 
binary coalescence events comes primarily from the 
discovery of a small number of galactic pulsars which 
appear to be in close binary systems with another neu- 
tron star. The first and most famous of these systems 
contains the Hulse- Taylor pulsar PSR B1316+16. It 
has demonstrated over many decades an orbital de- 
cay consistent with loss of energy to gravitational 

The number and inferred 



waves 



Weisberg et al. 2010 



lifetime of these systems can be used to obtain an es- 
timate of about one NS-NS merger event per Mpc 3 
per million years Abadie et al.|2010a , with up to two 



orders of magnitude uncertainty in rate due largely to 
our limited knowledge of the pulsar luminosity func- 
tion and limited statistics. The merger rate translates 
into an estimate of ^0.02 detectable NS-NS merger 
events per year for the initial LIGO- Virgo detector 
network, which operated between 2005-2010, and ~40 
per year for the advanced detectors once they reach 
design sensitivity. Although their gravitational radia- 
tion is stronger, the merger rates of NS-BH binaries is 
more uncertain as we have not observed any NS-BH 
binary systems, and have generally poor knowledge of 
the black hole mass distribution. 

Gamma-ray bursts (GRB) are flashes of gamma 
rays observed approximately once per day. Their 
isotropic distribution in the sky was the first evi- 
dence of an extra-galactic origin, and indicated that 
they were extremely energetic events. The duration 
of prompt gamma-ray emission shows a bi-modal dis- 
tribution which naturally groups GRBs into two cat- 
egories Kouveliotou et al. 1993] , Most long GRBs 
emit their prompt radiation over timescales >2 sec- 
onds, and as much as hundreds of seconds. They have 
been associated with young stellar populations and 
the collapse of rapidly rotating massive stars Hjorth 



and Bloom||2011| . Short GRBs (sGRB), with prompt 
emission typically less than Is and a generally harder 
spectrum, are found in both old and new stellar pop- 
ulations. Mergers of two neutron stars, or of neutron 
star/black hole systems, are thought thought to be a 
major contribution to the sGRB population Nakar 
2007]. It is this favored progenitor model which 
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makes short GRBs and associated afterglow emission 
a promising counterpart to gravitational-wave obser- 
vations. 

Since 2005, the Swift satellite has revolutionized our 
understanding of short GRBs by the rapid observation 
of x-ray afterglows, providing the first loca lization, 
host identification, and red-shift information Gehrels 
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et al. 2009 . The beaming angle for short GRBs is 



highly uncertain, although limited observations of jet 
breaks in some afterglows imply half-opening angles of 
3-14 degrees |Liang et al.|2008[|Fong et al.|2012 



The absence of an observable jet break sets a lower 
limit on the opening angle which is generally weak 
(due to limits in sensitivity), though in the case of 
GRB 050724A, late-time Chandra observations were 
able to constrain 9j > 25° Grupe et al.||2006 



The observed spatial density of sGRB's and lim- 
its on beaming angle result in a NS-NS merger event 
rate roughly consistent with that derived from galactic 
binary pulsar measurements. Although the beaming 
factor of ~ 0f /2 means we believe most merger events 
seen by the advanced GW detectors will not be accom- 
panied by a standard gamma-ray burst, this is some- 
what compensated by the fact that the ones that are 
beamed toward us have stronger gravitational-wave 
emission. Current estimates for coincident GW-sGRB 
observation for advanced LIGO- Virgo are a few per 



year assuming a NS-NS progenitor model Metzger 



and Berger 2012 Coward et al. 2012 . The rate in- 



creases by a factor of 8 if all observed short GRBs 
are instead due to NS-BH (10 M®) mergers which are 
detectable in gravitational-waves to about twice the 
distance. 

In addition to the jet-driven burst and afterglow, 
other EM emission associated with a compact merger 
can be a promising channel for GW-EM coincidence, 
particularly if the EM radiation is less-beamed or even 
isotropic. A few short GRBs (~10%) have shown 
clear evidence of high-energy flares which precede the 
prim ary burst by 1-10 seconds, and possibly up to 
100s Troja et aL||2010 . The precursors can be inter- 
preted as evidence of some activity during or before 
merger, such as the resonant shattering of NS crusts 
Tsang et al.||2012 , which could radiate isotropically. 



Thus it will be interesting to search for weak non- 
standard EM emission accompanying all nearby NS- 
NS mergers seen in GWs, while we expect only a small 
fraction to be oriented in our line-of-sight for a stan- 
dard jet-driven sGRB. 



The Gamma-ray Burst Monitor (GBM) |Meegan 



et al. 2009 aboard the Fermi spacecraft measures pho- 



ton rates from 8 keV-40 MeV. The instrument con- 
sists of 12 semi-directional Nal scintillation detectors 
and 2 BGO scintillation detectors which cover the en- 
tire sky not occluded by the Earth (about 65%). The 
lower-energy Nal detectors have an approximately 
cos 9 response relative to angle of incidence, and rela- 
tive rates across detectors are used to reconstruct the 
source location to a few degrees. The BGO detectors 
are much less directional, and are used to detect and 
resolve the higher energy spectrum above ~200 MeV. 

GBM produces on-board triggers for gamma-ray 
burst events by looking for multi-detector rate ex- 
cess over background across various energy bands and 
timescales. In the case of a trigger, individual pho- 



ton information is sent to the ground and the event 
is publicly reported. Those events which have been 
confirmed as GRBs have already been studied in co- 
incidence with LIGO- Virgo data Abbott et al.||2010 



Abadie et al. 2010b 2012a . So far, no gravitational 



wave counterparts to triggered GRBs have been iden- 
tified, which was not unexpected given the limited 
reach of the first generation instruments. 

In addition to the triggered events, survey data is 
available which records binned photon counts over all 
time. In this proposed offline analysis of short tran- 
sients, we consider the CTIME daily data, which con- 
tains counts binned at 0.256s over 8 energy channels 
for each detector. A new GBM data product (contin- 
uous TTE) was implemented in late 2012. It provides 
continuous data on individual photons with 2/is and 
128 energy channel resolution, which will further en- 
hance offline sensitivity to particularly short bursts. 

This paper discusses the possibility of using these 
offline GBM data products to follow-up gravitational- 
wave candidates in the advanced LIGO- Virgo era. In 
section [2] we describe the characteristics of a trigger 
provided by gravitational- wave data. Sections [3] and [4] 
develop a likelihood-ratio based procedure for analyz- 
ing available GBM offline data about the time and sky 
location provided by the gravitational waves. Finally, 
in section [5] we demonstrate the performance of the 
search algorithm on background times and a sample 
of known sGRB's. 



2. Gravitational-wave trigger 

Searches for NS-NS and NS-BH coalescence in 
gravitational-wave data typically use matched filter- 
ing of the model waveform Abadie et al. |2012b|c 



The gravitational waves from CBC are character- 
ized by a chirp of monotonically increasing frequency 
and amplitude as the binary inspirals under approx- 
imately adiabatic orbital decay. During this progres- 
sion, kilometer-scale interferometers are sensitive to 
the inspiral phase of a stellar-mass coalescence until 
just prior to merger, but before any tidal effects on a 
NS become important. 

The waveform in this regime is well-modeled us- 
ing post-Newtonian techniques, and can be consid- 
ered a standard candle encoding orbital parameters 
of the system. Mass, spin, and coalescence time are 
accessible with the GW projection onto a single de- 
tector, while distance, inclination, and sky location 
degeneracies can be disentangled to varying degrees 
of success using coherent data from multiple detec- 
tors, for example by using Bayesian inference Veitch 



et al. 2012 . A typical gravitational-wave detection 
in the early advanced detector era may have moder- 
ate signal-to- noise ~8 in two or more detectors, es- 
timated merger time to within milliseconds, and sky 
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simulated NS/NS merger localization 
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Figure 1 : Example sky localization of simulated NS-NS 
merger seen in modeled gravitational-wave data with 
noise similar in spectrum to that expected from the 
advanced LIGO- Virgo GW network. The signal-to-noise 
of the GW signal in each detector is 8.7 (HI) 7.2 (LI) 
and 3.1 (VI), and the sky location is determined through 
Bayesian inference using a nested sampling technique 
|Veitch et al.|2012 applied to data from all detectors. 
The elongated shape beyond 1<t is due to a small 
degenerate timing uncertainty in VI. 



localization to ^100 square d egrees dependin g on de- 
tector network configuration Fairhurst||2011 . Figure 
[T] shows a reconstructed 3-detector localization for a 
simulated signal near detection threshold. 



3. Coherent analysis of GBM data 

In this section, we develop a procedure to coher- 
ently search GBM detector data for modeled events. 
The basic idea is that by processing multiple detec- 
tor data coherently, we can obtain a greater sensi- 
tivity than when considering one detector at a time. 
Greater computational resources available offline (vs. 
on-board) also allow for more careful background es- 
timation to be done. For this analysis, we can relax 
to some extent the strict 2-detector coincidence re- 
quirement used to veto spurious events on-board as 
the gravitational-wave trigger means much less time 
and sky area is considered. 



3.1. GBM background estimation 

Each detector is subject to a substantial time- 
varying background from bright high-energy sources 
that come in and out of the wide field of view, as well 
as location-dependent particle and Earth atmospheric 
effects. This background must be estimated and sub- 
tracted out to look for any prompt excess. Methods in 



use include the local averaging of previous data done 
on-board, smooth spline-fits with a high-frequency 
cutoff, direct tracking and modeling of the dominant 
sources, and averaging rates from previous orbits with 



similar orientations Finger et al. 2009 


Wilson-Hodge 


et al.| 


2012 




Fitzpatrick et al. 


2012 . 


In this analy- 



sis where we are interested in the background esti- 
mate for a short foreground interval [— T/2, +T/2] 
where T ~ls, we estimate the background using a 
polynomial fit to local data from [— 10T, +10T] (min- 
imum ±5s), excluding time [— 3T/2, +5T/2] around 
the foreground interval to avoid bias from an on-source 
excess. An example of the foreground and background 
intervals about a strong prompt signal is shown in fig- 
ure [2] The polynomial degree is determined by the 
interval length to account for more complicated back- 
ground variability over longer intervals. It ranges from 
2 (minimum) to 1+0.5 log 2 T. The quality of the fit 
is determined using a x 2 statistic applied to the data, 
re-binned at T/4 seconds for T > Is. If the x 2 P er 
degree-of- freedom is over 2, or if any of the N indi- 
vidual data points has x 2 > 41niV, it is assumed that 
the polynomial could not adequately model the local 
background variability, and the fit is redone over the 
smaller interval [— 5T, +5T]. If the fit continues to fail 
the x 2 test & t the looser requirements of 3 and 6 In N 
respectively, the background estimate is marked un- 
reliable and not used in further computations for that 
particular foreground interval. 

High-energy cosmic rays striking a Nal crystal can 
result in long-lived phosphorescent light emission. 
The detector may interpret this is a rapid series of 
events, creating a short-lived jump in rates for one or 
multiple channels, and severely distorting the back- 
ground fit if not accounted for. They are identified 
with a simple procedure that compares the counts in 
each 0.256s bin against the mean rate estimated from 
four neighboring bins. If an excess is detected to 5a or 
more, that measurement and the immediately neigh- 
boring bin on either side is removed from the back- 
ground fit for all channels of the affected detector. 
A much stricter requirement is used when rejecting 
cosmic-ray events in the foreground interval. In that 
case, a 0.256s measurement must have signal-to-noise 
over 500 /y/T for it to be excised, which is much larger 
than expected for true GRBs. 

3.2. Likelihood-ratio statistic 

A likelihood ratio combines information about 
sources and noise into a single variable. It is defined 
as the probability of measuring the observed data, d, 
in the presence of a particular true signal Hi (source 
amplitude s > 0) divided by the probability of mea- 
suring the observed data in noise alone H (s = 0). 



A(d) 



P(d\H ) 



(1) 
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Figure 2: Detector and channel-dependent polynomial 
background fitting applied to GRB 090510A. Four of the 
8 x 14 channel and detector combinations are shown. 
The green represents the highest likelihood foreground 
interval, and the background contribution is estimated 
using a polynomial fit to the region (cropped) shown in 
blue. Red dots mark the average flux in the foreground 
interval. 



free sky-location parameters in Vj and the single am- 
plitude parameter s. The shape of the likelihood func- 
tion over source amplitude s is the product of (almost) 
Gaussians centered about the best-estimate of s in 
each measurement, and the total likelihood is scaled 
by the sum of the squared signal-to-noise of all mea- 
surements. 

The dependence of response factors r$ on sky loca- 
tion is complicated, so the likelihood ratio is calcu- 
lated over a sample grid of all possible locations. As- 
suming a single location, the remaining free param- 
eter is the source amplitude s. The variance in the 
background-subtracted detector data includes both 
background and source contributions, 



+ ris - 



2 2 
0"„.S 



( s >0) 



(5) 



with a 2 , representing Gaussian-modeled systematic 
uncertainty in the instrumental response. Source 
terms are only included for physical s > else their 
contribution is zero. The background contributes 
Poisson error, as well as any systematic variance a 2 . 
from poor background fitting which is also assumed 
to be Gaussian, 



When signal parameters such as light-curve, spec- 
trum, amplitude s and sky-location a, S are unknown, 
one can either marginalize over the unknown param- 
eters, or take the maximum likelihood over the range 
to obtain best-fit values. 

For i measurements of sufficiently binned, uncorre- 
cted Gaussian data, 



(di - 7-jS) 2 



(2) 



(3) 



where we have used di = di — (rij) to represent the 
background-subtracted measurements (e.g. red dots 
minus blue curve in figure |2J), o~ ni and for the stan- 
dard deviation of the background and expected data 
(background+signal) , n for the location/spectrum- 
dependent instrumental response, and s, a single in- 
trinsic source amplitude scaling factor at the Earth. 
Maximizing the likelihood ratio is the same as maxi- 
mizing the log-likelihood ratio C = In A, 



In 



°~di 



d~i 



(d, -r lS ) 2 
2al 



(4) 



When the first two terms are fixed, maximizing the 
log-likelihood ratio is equivalent to minimizing the the 
third term which we recognize as a x 2 fit to the two 



(0) 



We find Sbcst which maximizes C by setting the deriva- 
tive dC/ds to zero. If a^. = a n ., which happens when 
rs of, V2£max reduces to a coherent SNR (sum 
data using weights ri/a 2 .). If can be assumed 
constant, the solution for Sbest is found analytically, 



^best 



Ei Tidily 



(7) 



which is just an appropriately inverse-noise weighted 
sum of the individual estimates of s from each mea- 
surement. Although depends on s, we can make 
the practical approximation er^ w max((rij), d{) + o~\. 
as an initial guess. 

To find the true maximum, we solve for dC/ds = 
using Newton's method, beginning with the value of 
the first derivative at our initial guess sq, and using the 
analytic second derivative to refine the measurement: 
Sbest ~ so — (dC/ ds) j (d 2 C j ds 2 ) . This calculation is 
both fast and easily vectorized. One initial guess using 
equation [7] followed by a couple iterations of Newton's 
method generally provides an excellent approximation 
for Sbcst- 

So far we have been calculating the probability of 
a signal assuming a specific source amplitude s. To 
consider all possible source amplitudes we need to in- 
tegrate the likelihood P(d\s) (equation [5]) over a prior 
on s, 



P(d) = / P(d\s)P(s)ds 



(8) 
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For a given set of detector data d, the likelihood P{d\s) 
over s is almost the product of individual Gaussian 
distributions (not quite Gaussian because o~d depends 
on s). The product of Gaussian distributions with 
mean values /ii and standard deviations cr, is itself 
Gaussian with mean and variance, 



Mpiod 



E mtf 



i 



E i/of 



(9) 



In this case, /Uj = dj/ri and <7j = ov^/r^. The mean 
value /iprod is the same as our initial guess for s at 
maximum likelihood (equation [7]) , but we can assume 
to have a more accurate maximum-likelihood location 
Sbest from the numerical procedure outlined above. 
The estimate for the variance of C over s is, 



4 



i 



E^ 2 



o-\. evaluated at 

^bcst 

(10) 



For a flat prior P(s) = 1, we can integrate equa- 
tion [8] by simply considering the area of a Gaussian 
with peak value P(d|sbest) and variance crj. Other 
choices for a prior can be represented by a power-law 



distribution P(s) oc 



-P 



A spatially homogeneous 



population, suitable for nearby sources, would assume 
ft = 5/2, while the empirical distribution for observed 
GRB amplitudes follows a power-law decay closer to 
j3 1.8, reflecting cosmological effects. If instead we 
arc looking for signals from a particular host galaxy 
at known location and distance, we want to use an 
intrinsic source luminosity distribution. A convenient 
option is to use a scale-free prior with fixed j3 = 1, 
so that our choice of form for an amplitude prior at 
the Earth does not translate back into a luminosity 
distribution that varies with distance. 

One difficulty with any power-law prior is that it 
diverges for s — > 0. In a reasonable scenario of SNR 
>a few, the integral over s will consist of two distinct 
contributions. The first is a Gaussian component with 
t 1 = Sbcst and variance cj, scaled by the prior P(sbcst)- 
The second component is an infinite contribution from 
the divergent prior at s — with little contribution 
afterward due to suppression from the Gaussian tail. 
The infinite contribution represents the certainty of a 
signal of arbitrarily small amplitude to be present in 
the data, regardless of the ability of the data to resolve 
it. For moderate SNR, it's very easy to isolate only 
the Gaussian contribution by placing a small cut on 
amplitude, truncating the likelihood for s < s cut . Un- 
der the approximation that P(s) varies slowly over the 
width of the Gaussian, the log-likelihood marginalized 
over all s becomes, 



£(d) 



f In Sbest + In o- c +C{d\ s best ) (11) 



up to common additive constants that do not depend 
on the data d. 



The approximation for the marginalized likelihood 
becomes problematic for small Sbcst where the as- 
sumption that P(s) is constant over the range of the 
Gaussian breaks down, and the Gaussian distribution 
can no longer be isolated from the divergence at small 
s. We enforce a finite and well-behaved prior by mul- 
tiplying by a prefactor, 



P(s) = [l - e-^l^f 



(12) 



so that P(s) reaches a maximum constant value of 
(l/7°£) /3 f° r small s. The tunable parameter 7 sets 
the number of standard deviations at which the prior 
begins to plateau, and we use 7 = 2.5. This allows 
us to use the approximation that P(s) is reasonably 
constant over a range of o£ for any s > 0. The only 
remaining correction is to account for clipping of the 
Gaussian for non-physical s < 0, which can be rep- 
resented by the error function. The final approxima- 
tion for the amplitude-marginalized log-likelihood be- 
comes, 




1 + Erf 

I _ e -(sbc 5 t/7^c)' 3 

-/3 In (ja c ) 



^best 



+ £(d|s be st) 



V2ac 

- P In Sbcst Sbcst > 

Sbcst < 



(13) 



which contains factors form the Gaussian width, frac- 
tional overlap with s > 0, maximum likelihood at 
Sbest, and scaling from P(s) respectively. Finally we 
are free to calibrate the log-likelihood by subtracting 
the expected C(d) calculated for no signal at a ref- 
erence sensitivity: £ rc f = — ,3 In 7 + (1 — 0) lncr£ jrc f . 
o£ represents the source amplitude required for a 
lcr excess in the combined data, and is around 0.05 
photons/s/cm 2 x (T/ls)" 1 / 2 [50-300 keV] for typi- 
cal source spectra and reference background levels. 
Figure [3] shows the coherent signal-to-noise expected 
from all detectors for a 0.512s-long event with nor- 
mal GRB spectrum and constant amplitude of 1.0 
photons/s/cm 2 , and compares it to the SNR expected 
from the most favorably-oriented individual detector 
alone in the 50-300 keV band. 



4. Performing the follow-up 

For a single gravitational-wave trigger, we search 
a standard sGRB accretion timescale of [0, 5s] rel- 
ative to the time of coalescence for prompt flux ex- 
cess in GBM between 0.256 and 2s long (the lower 
limit of 0.256s set by the CTIME archival resolution 
does not apply to continuous TTE data). We also 
search a prior interval [-30s, 0s] for possible precur- 
sor bursts between 0.256s and 2s. Finally we may 
include bursts between 2s to 32s in an extended inter- 
val [-30s, 300s] to search for possible longer-duration 
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by beginning with the highest-ranked event and re- 
moving from the list any lower-ranked events with 
foreground windows that overlap it; then taking the 
next surviving highest ranked event, etc., until the list 
is exhausted. 



~o 5. Test on Swift short GRBs 



Figure 3: Sigrial-to-noise expected from GRBs with a 
normal spectrum, and coming from various positions 
across the GBM field-of-view. The hypothetical signal 
lasts 0.512s and is normalized to 1.0 photons/s/cm 2 in 
the 50-300 keV band, while background rates and Earth 
position are selected from 10 seconds prior to GRB 
090305A. The instrumental response includes 
contributions from atmospheric scattering in the 50-300 
keV band. The first map shows the signal-to-noise from 
the single best detector over 50-300 keV. The second plot 
shows the SNR expected from a coherent analysis of the 
full CTIME data. The coherent statistic must be 
calculated for all directions, which means it is subject to 
a trials factor equal to the number of independent sky 
locations. 



emission. While emission outside of the standard ac- 
cretion timescale can be considered speculative, it's 
worth looking for given that any events detectable in 
gravitational-waves will be closer than sGRB's with 
known red-shift to date, making them good candidates 
to search for weak exotic and possibly less-beamed 
emission. To appropriately tile the search in time and 
duration T, we use rectangular windows with T spaced 
by powers of two (0.256s, 0.512s, Is, etc.). Their cen- 
tral times are sampled along the search interval in 
units of T/4 to provide an even mismatch in signal- 
to-noise across search windows. 

The likelihood-ratio statistic described in the pre- 
vious section is calculated for each foreground inter- 
val based on the background-subtracted flux measure- 
ments di in each of the 8 x 14 channel and detec- 
tor combinations. In addition, the instrumental re- 
sponse ti (and thus the likelihood-ratio) also depend 
on source sky location and spectrum. To efficiently 
calculate the response over a large area of sky, we use 
precomputed all-sky response loo k-up tables originally 
generated for offline localization Paciesas et al.|2012~ 



The likelihood ratio as a function of sky position pro- 
vides a probability distribution over the sky of a GBM 
signal (as in figure [5]) . This can be coincided with the 
GW-derived skymap (figure [I]) by direct multiplica- 
tion, in effect using the GW skymap as a prior. 

Finally, we marginalize over sky location and rep- 
resentative source spectra to get a ranking of events 
characterized only by their foreground time intervals. 
A unique list of non-overlapping events is constructed 



We can test the offline analysis on short GRBs trig- 
gered by Swift and observable by Fermi-GBM. The 
Swift GRBs are particularly useful as their accurate 
localizations resolve systematic errors in the GBM 
model response. For this test, we calculate the all- 
sky marginalized likelihood ratio for GBM data local 
to Swift sGRB measurements (figure [4]) . Here we as- 
sume no previous information about source location. 
A background rate distribution is obtained by running 
the same search over nearby time. 

About 50% of Swift sGRB's (T90 < 2s) are within 
GBM's field-of-view (the 65% of the sky not occulted 
by the Earth) and occur during an interval of time 
when Fermi is outside of the South- Atlantic Anomaly 
(SAA) and operational. Most of these observable 
Swift sGRB's also triggered GBM on-board. GRB 
081024A triggered on-board, but was too close to the 
interruption of data taking during passage through the 
SAA for our two-sided background estimation. GRB 
090305A and 120403A did not trigger on-board, and 
do not show compelling evidence for a signal in the 
offline data. GRB 090305A and 120403A also did not 
trigger on-board, but do show clear evidence of a sig- 
nal present in the offline analysis (figure [5| , though 
with relatively poor statistics. The remaining Swift 
short-GRBs from the observable sample both trig- 
gered on-board and are clearly identified offline. 



6. Conclusion 

Direct detection of gravitational-waves from the 
merger of NS-NS or NS-BH binary systems is ex- 
pected to occur within the next few years as advanced 
2nd generation ground-based gravitational-wave de- 
tectors come online. Fermi GBM provides a unique 
and promising opportunity to observe EM counter- 
parts due to its large sky coverage, and the anticipated 
association between NS-NS mergers and sGRB's. To 
maximize the use of information from both GWs and 
GBM, we have outlined a joint search strategy trig- 
gered by observations of GW signals from coalescing 
binary systems in which a small amount of local GBM 
data is scanned using a likelihood-ratio based analy- 
sis applied to the full instrument data. The sensitivity 
of the method to short GRBs will benefit in the fu- 
ture from continuous GBM offline data products with 
higher time resolution and improved all-sky response 
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Figure 4: The maximum likelihood coherent flux excess 
found by searching time around Swift short-GRBs within 
the GBM field of view using foreground intervals between 
0.256 and 2s. GRBs are placed at 0.01 Hz in anticipation 
of a GW-GBM coincidence search window of <100s. An 
estimated cumulative background rate distribution is 
found by taking the non-overlapping all-sky events from 
a continuous scan of approximately one-day in total of 
nearby data. The dotted background distribution shows 
the effect of a simple cut to remove particle events (the 
cut does not remove any GRBs). Of the four 
un-triggered events, GBM excesses corresponding to 
GRB 090815C and 110112A are buried within the 
background distribution, and are not clearly associable 
with their GRBs. The ability of the method to identify 
the next two weak, un-triggered events GRBs 090305A 
and 120403 A (see also figure [5]l depends on ongoing 
strategies to further reject non-Gaussian outliers. The 
remaining Swift sGRB's in the sample have a 
corresponding GBM on-board trigger. 



models. We anticipate such an effort will allow sen- 
sitive follow-up of NS-NS and NS-BH mergers, most 
of which are expected to not be accompanied by trig- 
gered GRB observation. 



090305A 
® Earth 



SNR: 8.9 
Chi-sq: 1.2 
"Spectrum: n 



Fermi sky at 120403 01:05:22.288 +2. 048s 

50 100 150 200 250 300 350 „ 



» 120403A 
e Earth - 20 




SNR: 9.5 
Chi-sq: 1.2 
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Figure 5: Swift short GRBs (T90 < 2s) 090305A and 
120403 A were seen by Swift-BAT but did not trigger 
GBM on-board. Low photon statistics for these weak 
events result in large sky location uncertainties shown by 
the 1, 2, and 3cr confidence regions. The quoted 
signal-to-noise ratio is calculated at the best-fit location 
using the full spectral data from all Nal and BGO 
detectors. The spectrum represents one of three 
representative spectral models (hard, normal, soft) for 
the model response that gives maximum likelihood, while 
the skymap is marginalized over all spectral models. 
Arrows represent movement of the Fermi spacecraft over 
±2 minutes. 
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